Magnetization plateaus in spin chains: "Haldane gap" for half-integer spins 
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We discuss zero-temperature quantum spin chains in a uni- 
form magnetic field, with axial symmetry. For integer or half- 
integer spin, S, the magnetization curve can have plateaus 
and we argue that the magnetization per site m is topologi- 
cally quantized as n{S — m) = integer at the plateaus, where 
n is the period of the groundstate. We also discuss conditions 
for the presence of the plateau at those quantized values. For 
S = 3/2 and m = 1/2, we study several models and find two 
distinct types of massive phases at the plateau. One of them 
is argued to be a "Haldane gap phase" for half-integer S. 

PACS numbers:75.10.Jm 



One-dimensional antiferromagnets are expected not to 
have long-range magnetic order in general. It was argued 
by Haldane jl), in 1983, that for integer, but not half- 
integer spin, 5*, there is a gap to the excited states. In 
the presence of a magnetic field, the 5=1/2 Heisenberg 
antiferromagnetic (AF) chain remains gapless from zero 
field up to the saturation field, where the groundstate is 
fully polarized ||]. For integer 5, the gap persists up to a 
critical field, equal to the gap, where bose condensation 
of magnons occurs ||^ . The S — 1 Heisenberg AF chain is 
known to be gapless from the critical field up to the satu- 
ration field Recently Hida observed that an S — 1/2 
antiferromagnetic chain with period 3 exchange coupling 
shows a plateau in the magnetization curve at magnetiza- 
tion per site m = 1/6 (1/3 of the full magnetization) 
Related works on bond-alternating chains have also been 
reported |^-^ including experimental observation pd| ]. 
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FIG. 1. Schematic diagram of possible magnetization curve 
of an S = 2 chain. There are plateaus at special values of 
magnetization m. The plateau with a fractional S — m = 2/3 
must accompany a period 3 groundstate. (See text for details.) 



In this letter, we consider the zero-temperature be- 
haviour of general quantum spin chains, including chains 
with periodic structures, in a uniform magnetic field 
pointing along the direction of the axial symmetry (z- 
axis). (i.e. the total is conserved.) We argue that, 
in quantum spin chains, there is a phenomenon which is 
strikingly analogous to the Quantum Hall Effect - topo- 
logical quantization of a physical quantity under a chang- 
ing magnetic field [|l^.(See Fig. [^.) We first consider an 
extension of the Lieb-Schultz-Mattis (LSM) theorem jl^] 
to the case with an applied field. This indicates that 
translationally invariant spin chains in an applied field 
can be gapful without breaking translation symmetry, 
only when the magnetization per spin, m, obeys S — m — 
integer. We expect such gapped phases to correspond 
to plateaus at these quantized values of m. "Fractional 
quantization" can also occur, if accompanied by (explicit 
or spontaneous) breaking of the translational symmetry. 
The generalized LSM theorem does not prove the pres- 
ence of the plateau, however. Thus we construct a cor- 
responding argument using Abelian bosonization, which 
is in complete agreement with the generalized LSM the- 
orem, and also gives a condition for the presence of the 
plateau. As simplest examples, we study translationally 
invariant 5 = 3/2 chains at m = 1/2. We present nu- 
merical diagonalization and Density Matrix Renormal- 
ization Group (DMRG) ||T^ calculations, which demon- 
strate the existence of the two distinct types of gapped 
phases for generalized models. They are related to the 
5 = 1 large- phase and the 5=1 Haldane phase, re- 
spectively. On the other hand, our study shows that the 
standard 5 = 3/2 Heisenberg model is gapless with no 
plateau at to = 1/2. We only give a brief summary of our 
numerical work in this letter; details and further results, 
including the effect of the axial symmetry breaking, will 
be presented in a longer paper JTst . 

The LSM theorem |]l^ proves the existence of at least 
one low-energy, 0{1/L) excited state for even length L 
half-integer 5 AF chains with periodic boundary condi- 
tions and general, translationally invariant Hamiltonians. 
It is expected that, this implies either gapless excitations 
or spontaneously broken translational symmetry in the 
L oo limit. The failure of this proof for the integer-5 
case is necessary for the existence of the Haldane phase 
with no broken translational symmetry and a gap. We 
observe that the original version of this proof also works 
in a magnetic field except for integer S ~ m. Thus only 
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in this case is a massive phase without spontaneously 
broken translational symmetry possible. The proof con- 
sists of making a slow rotation on the groundstate, li/i), 
assumed to be unique, and observing that the resulting 
low-energy state is orthogonal to the groundstate. The 
rotation operator is U = exp [— i^^^^(27rj/L)5'|]. For 
any Hamiltonian H, including a magnetic field term, with 
short-range interactions which is invariant under rotation 
about z-axis, it can be shown that 
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This implies the existence of an excited state with ex- 
citation energy of 0{1/L), if we can show that U\ip) is 
orthogonal to \^). To this end, we use the invariance of 
H under translation by one site: T. This operation maps 
U into: 



U 



(2) 



Namely, the operation of U changes the eigenvalue of T 



Thus 



by a factor e^^'^^^""), where m = J2j=iS-/L. 
Ultp) must be orthogonal to {ip) except when (S — m) is 
an integer. We note that this is consistent with previous 
results for translationally invariant S' = 1/2 and 1 AF 
chains, where no gap is found at partial magnetization. 
However for higher spin, gapped phases at partial magne- 
tization are possible without breaking the translational 
symmetry, when — to is an integer. 

When S* — TO is not an integer, there is a low- lying state 
with energy of 0(1/ L). This means either a massless 
phase with a continuum of low-energy states or sponta- 
neous symmetry breaking in the thermodynamic limit. 
Following the above proof, when S — m = p/q where p 
and q are coprimes, U''\ip) for k = 0,1, . . . ,q — 1 have 
different eigenvalues of T. Thus these q states have low 
energy of 0{1/L). If these are related to a spontaneous 
breaking of the symmetry, the ground states in the ther- 
modynamic limit should be g-fold degenerate. Since they 
have q different eigenvalues of T, they can be related to 
a spontaneous breaking of the translation symmetry to 
period of q sites in the thermodynamic limit. It is natural 
to expect a gap and plateau in this case. As in the case of 
Quantum Hall Effect, "fractional quantization" is there- 
fore possible, accompanying the spontaneous breaking of 
the translation symmetry in the present case. (See Fig. ^) 
We may compare this to a hidden symmetry breaking in 
Fractional Quantum Hall Effect |16|. 

Possibly there is a low-energy state other than con- 
structed as above. In such a case, we can again con- 
struct a set of q low-energy states by operation of U. 
Thus in general the period of ground state can be an in- 
tegral multiple of q. A simple example of such a case 
is known for 5 = 1 and m = 0. While a gap with- 
out spontaneous symmetry breaking is possible as in the 
Haldane phase, spontaneous dimerization occurs in the 



bilinear-biquadratic Heisenberg model for a range of pa- 
rameters PH . 

Our generalization of LSM theorem is easily extended 
to Hamiltonian with spatial structures: bond-alternating 
chains [^,D, spin-alternating chains px| ], spin ladders, 
etc. For example, Hida's model fq] is only invari- 
ant under a three-site translation T'^; a massive phase 
without spontaneous symmetry breaking is possible for 
3(5 — to) = integer. Thus a quantized plateau is possible 
at TO = 1/6 as he observed. In general, the quantiza- 
tion condition is given by Su — to^ = integer, where Su 
and TTiu are respectively the sum of 5 and to over all 
sites in the unit period of the groundstate. The period 
of the groundstate is determined by the explicit spatial 
structure of the Hamiltonian, and also by spontaneous 
symmetry breakings. 

The low-energy state Ultp) appearing in the LSM the- 
orem has the same total magnetization as in the ground- 
state. It does not directly contradict the existence of a 
plateau, which is determined by the gap to states with 
other total magnetizations. However, we expect that, 
in general a gapless phase has low-energy states in both 
fixed and different magnetization sector, as can be seen in 
the following Abelian bosonization approach. Schulz |l^] 
explained the difference between integer and half-integer 
spin by Abelian bosonization. We show that his result 
can be understood more simply as a consequence of sym- 
metries. At the same time, we generalize the discussion 
to the case with a non-vanishing magnetization. Fol- 
lowing Schulz p7[ |, we start from Abelian bosonization 
of 25 spin- 1/2 chains and then couple them to form a 
spin-5 chain. Firstly, each spin-1/2 chain is fermion- 
ized by Jordan- Wigner transformation. The z compo- 
nent of each spin-1/2 is related to the fermion number 
as (T^ = 1 — 2ip'l^tpn- Then the low-energy excitations are 
treated by continuous fermion fields. Let us denote the 
lattice spacing as a and the spatial location x = na. The 
continuous fermion fields ippi. and ipL are defined by 
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where j = 1, • • • , 25 is the "flavor" index to distinguish 
25 spin-l/2's. They are bosonized in a standard way: 
-0;^ — e^'^R^^ and ■(/'i = e~*'''i/"'^, where i^-^ and tp-'j^ are 
chiral bosons and R is the compactification radius of the 
boson. R will be renormalized by interactions [Q, and 
will eventually depend on the model and on the magne- 
tization TO. (For an isotropic model, R is fixed by the 
symmetry at m = 0, but the magnetic field breaks the 
symmetry and thus R will depend on to.) We define 
the non-chiral bosonic field (^^ = -t- LfPj^ and its dual 

Interactions among bosonic fields are also generated 
during the mapping from the original spin problem. In 
general, we expect any interaction would be generated if 
not forbidden by a symmetry. Thus we analyze symme- 
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tries of the system, following the treatment of spin-1/2 
chains in Ref. The original problem has a C/(l) 

symmetry: rotational invariance about z axis. Rotation 
of each spin-1/2 is given by the phase transformation 
■0^ ^ — > e*^V'i B. of the corresponding fermion. In bosonic 
language, this corresponds to a shift of the dual field 
Cpi (pj + const.. Since we have coupled 25* spin-1/2 
chains into a spin-S* chain, only the simultaneous rota- 
tion of 2S spin-l/2's is a symmetry of the system. If 
we define a new bosonic field (f) — if^ (and similarly 

for (f^), the U{1) symmetry is written as ^ —*■(/> -I- const. 
Thus all the interactions of the form e^^niriiiii!. prohib- 
ited by the symmetry. The remaining 25" — 1 fields, which 
are defined by linear combinations of original ip^ fields, 
are not protected by the symmetry. Thus all fields ex- 
cept (j) are expected to become massive by interactions, as 
Schulz observed by an explicit calculation. The remain- 
ing (j) field, is also subject to e^™*^/^ type interactions. 
Let us consider another symmetry of the system: one-site 
translation. By definition (^, it actually corresponds to a 
transformation of the continuum field ^/ijj — > e^'^'^'^i/'j^ ^^'^ 
4>l^ e~^^'"^4'L- Again, only the simultaneous transla- 
tion of all flavors is a symmetry of the system. Thus the 
one-site translation T is written as -I- 4:S{kpa)TTR, 

in the bosonic language. 

Since all of the 25* flavors are equivalent, the mag- 
netization should be equally distributed among them. 
Thus the fermi momentum kp is determined as kpa = 
{S — m)n/{2S) . As a consequence, the one-site transla- 
tion T is given by 

(f) (I) + 2{S - m)TrR. (4) 

Thus the leading operator cos {4>/R) is permitted only if 
S' — m is an integer. For m satisfying the quantization 
condition, the leading operator cos {4>/R) should be rel- 
evant in order to produce a gap. Thus R must be larger 
than Rc = I/VStt for the presence of the plateau. If 
S — m — p/q where p and q are coprimes, the operator 
cos {q<p/ R) is permitted. It can be relevant if i? > q/^/Sn 
(this is a severe condition for a large q); if it is, a ground- 
state in the thermodynamic limit corresponds to a po- 
tential minimum of cos {q(j)/ R). There are q such ground- 
states and they are mapped to each other by applying the 
translation operator T*' {k < q). Thus the ground states 
have spontaneously q-fold broken translation symmetry. 
These results are in agreement with the generalized LSM 
theorem, and also give conditions for a finite plateau at 
the quantized values. 

Our bosonization argument is also readily generalized 
to models with spatial structures. Our picture is con- 
sistent with Okamoto's analysis [|j of Hida's plateau 
For S = 1 AF chains, Tonegawa et al. obtained an 
m = 1/2 plateau as soon as they introduced a small 
bond-alternation. In our approach, the leading opera- 
tor is expected to appear as soon as the translational 



FIG. 2. The partially magnetized VBS state for 5 = 3/2. 
A solid line denotes a valence bond (singlet formed from two 
spin-l/2's). An up arrow denotes a spin-1/2 with 5*^ = 1/2. A 
dashed circle represents the symmetrization of spin-1/2 vari- 
ables at each site. 



symmetry is broken. Thus we expect a plateau for any 
finite amount of bond-alternation, if the radius exceeds 
the critical value, in agreement with Ref. This is also 
in agreement with an explicit bosonization calculation by 
Totsuka |9) for S' = 1 bond-alternating chains. 

Now let us discuss some examples of translationally 
invariant 5 = 3/2 chains. It is interesting both from 
an experimental and conceptual point of view to add an 
easy-plane crystal field term: 

H = J2Sj- Sj+i + D{S;' f - hS]. (5) 
i 

Clearly, if >> 1, all the spins are first fixed to S*^ = 1/2 
with increasing field before any of the spins go into the 
= 3/2 state, corresponding to a gapped m = 1/2 
plateau. The presence of finite gap and plateau is proved 
rigorously for a sufRciently large but finite D [|l^, by 
applying the general theorem in Ref. [ pO[ . This situation 
is reminiscent of that which occurs in the large- Z? phase 
in a zero field 5 = 1 chain. Numerically, we found a finite 
m = 1/2 plateau at least for D >2. 

Another kind of trial groundstate for an S = 3/2 
chain, corresponding to an m — 1/2 plateau is shown 
in Fig. H in the valence bond notation |2^ . Regarding 
each 5 = 3/2 operator as being a symmetrized prod- 
uct of three S — 1/2's, one S — 1/2 is polarized by 
the applied field at each site while the other two form 
a valence-bond-solid (VBS) groundstate as occurs for an 
5 = 1 chain in zero field. In Ref. ||2^, this partially mag- 
netized VBS state was proposed and the relevance to the 
magnetization process was suggested. A generalization 
of this kind of VBS-type state and further analysis were 
later done in Ref. |23). (See also fl^.) Clearly this sort 
of VBS-type state exists for all 5 and m such that S — m 
is an integer. 

We can construct a model to realize the 5 = 3/2 VBS- 
type state in Fig. || as a ground state: 

H = Y^ ^i^'"^^^ + <^Sj ■ Sj+i - hS'j, (6) 
j 

where Pg"''"''''^'' is the projection operator onto the space 
with total spin 3 for sites j and j + 1. At a = 0, any state 
constructed with one valence bond between neighboring 
sites is a groundstate. The ground state is thus infinitely 
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degenerate due to the "free" spin-1/2 at each site. Ap- 
plying an infinitesimal magnetic field, the degeneracy is 
lifted and the ground state is the above mentioned VBS- 
type state (Fig. ||). Thus the model with a = has 
an TO = 1/2 plateau starting from zero magnetic field. 
Turning on the Heisenberg term, a, the degeneracy at 
h = is lifted and a finite magnetic field is required to 
reach to — 1/2. For small value of a, however, we might 
still expect a finite to = 1/2 plateau. We studied this 
model with periodic boundary conditions by numerical 
diagonalization for up to 12 sites and found the to = 1/2 
plateau exists at least for a < 0.06. In contrast to the 
plateau at large positive D, which is related to the large- 
D phase in 5* = 1 chains, it is natural to relate this state 
to the 5 = 1 Haldane phase. 

For S* = 1, the Haldane phase is known to be dis- 
tinct from the large-D phase; these two massive phases 
are separated by a critical point Dc, where the gap van- 
ishes [ p7p4|| . The Haldane phase is characterized by 
the existence of a topological long-range order and 
gapless edge excitations in the open boundary condi- 
tions |2^]. These are understood as consequences of a 
hidden symmetry breaking pc| ]. One might suspect that 
the two types of S* = 3/2 massive phases at the to = 1 /2 
plateaus discussed above, correspond to distinct phases. 

If they are distinct, there should be a phase transition 
between them. In terms of Abelian bosonization, this 
phase transition may be understood as the vanishing of 
the coefficient of the allowed relevant operator cos {(j)/R), 
as in the case of 5 = 1 ||l^. We numerically measured 
the gap (width of the plateau) for the model: 

H = J2 o^S, ■ S,+, + DiS^f + Pi-'-'^+i) - (7) 

which interpolates between (||) and (^. For a = 0.03 
(fixed), we find the plateau vanishes at ~ 4.5, separat- 
ing the "Haldane gap" type plateau and the "large-D" 
type plateau. Moreover, we compared the spectrum at 
a = 0.03 and D = Q between open and periodic bound- 
ary conditions, and found evidence for edge states. In 
the large-D region, there are no such edge states. These 
indicate that the "Haldane phase" at to = 1/2 plateau, 
which accompanies the edge states, is distinct from the 
"large- 1? phase" . 

We also numerically examined the standard S* = 3/2 
Heisenberg AF chain with open boundary conditions, by 
DMRG up to 100 sites. We did not find an to = 1/2 
plateau in this case, in agreement with Refs. We 
emphasize that the absence is not a priori obvious. As 
we have shown, in terms of the free boson theory, the 
plateau would be present if the compactification radius 
R is greater than the critical value Rc and the coef- 
ficient of the most relevant operator cos {(f>/R) is non- 
vanishing. We have determined the compactification ra- 
dius from the spectrum for the open boundary condition 



obtained by DMRG, as i? = 0.95i?c < Rc- We note 
that R is rather close to the critical value, and possi- 
bly we can realize a plateau by a small perturbation of 
the standard Heisenberg Hamiltonian ||l^. On the other 
hand, while the radius is not completely well-defined at 
the massive D = 2 plateau, a similar analysis gives the 
estimate R ^ 1.2Rc > Rc- This result is consistent with 
the presence of the plateau. 

The plateaus that we have found are closely related |l^] 
to Mott insulating (or charge density wave) phases in 
models of interacting fermions or bosons |2^j2^. Simi- 
larly to those cases [^,^, we have found that the sin- 
gular part of the magnetization curve near a plateau is 
proportional to y^\h — hc\ where he is the critical field at 
(either) edge of a plateau, at least for examples we have 
studied. Our approach will also give new insights into 
models of interacting fermions or bosons [p^ . 
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Note added (January 1997): After the submission 
of the present letter, we received a preprint by Totsuka, 
which is a substantial enhancement of his presentation |^ 
at Japanese Physical Society meeting Fall 1996, and con- 
tains some of our general argument using bosonization. 
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